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Abstract 

We supplement (and critically overview) the existing extensive analy- 
sis of antiferromagnetic solution for the Hubbard model with a detailed 
discussion of two specific features, namely (i) the evolution of the mag- 
netic (Slater) gap (here renormalized by the electronic correlations) into 
the Mott-Hubbard or atomic gap, and (ii) a rather weak renormalization 
of the effective mass by the correlations in the half-filled-band case, which 
contrasts with that for the paramagnetic case. The mass remains strongly 
enhanced in the non-half-filled-band case. We also stress the difference 
between magnetic and non-magnetic contributions to the gap. These re- 
sults are discussed within the slave boson approach in the saddle-point 
approximation, in which there appears a non-linear staggered molecular 
field due to the electronic correlations that leads to the appearance of 
the magnetic gap. They reproduce correctly the ground-state energy in 
the limit of strong correlations. A brief comparison with the solution in 
the limit of infinite dimensions and the corresponding situation in the 
doubly-degenerate-band case with one electron per atom is also made. 
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1 Introduction 

Antiferromagnetism (AF) appears in the Hubbard and related models for an 
arbitrary interaction strength U provided we are close to the half-fillcd-band 
situation (n — * 1 in the orbitally nondegenerate case) [T|. This is easy to un- 
derstand qualitatively, since the intraatomic Hubbard interaction Uj^i ( n «T n u) 
is diminished by keeping apart the electrons with the opposite spins. In the 
strong-correlation limit [2J this interaction leads to an antiferromagnetic kinetic 
exchange for an arbitrary band filling n. At the same time, the band energy is 
not increased because the concomitant nesting condition for the quasiparticle 
states achievable for bipartite lattices only leads to the energy decrease of the 
occupied states, even when going beyond the Hartree-Fock picture. In effect, 
the regime of the band filling n, for which the AF state is stable at given U has 
been determined for variety of theoretical approaches PP - 0]. The reliability 
of the results for the half-filled case is not in question as long as they reduce to 
those in the Hartree-Fock and to the mean-field (Heisenberg) approximations 
in the weak- and strong-correlation limits, respectively, as we will discuss in the 
following. The theoretical results are in accord with the fact that all known 
Mott insulators with the half-filled band configuration are also antiferromag- 
netic insulators. The main purpose of this paper is to overview the situation 
by discussing the crossover behavior from the Hartree-Fock to the mean-field 
approximation for the Heisenberg antiferromagnet. Similar, though not equiv- 
alent results can be reached within the dynamic mean-field theory (DMFT, see 
below). In the orbitally degenerate situation (with the degeneracy d > 1) the 
same type of magnetic(Slater) gap is generated by an alternant orbital ordering 
PI, as is also discussed at the end of the paper. 

Explicitly, we concentrate our attention to two specific features of quasiparti- 
cle states not elaborated in detail so far, namely, (i) an evolution of the magnetic 
gap (renormalized by the electronic correlations) into the Mott-Hubbard gap, 
and (ii) a rather weak renormalization of the effective mass for the half-filled- 
band case, which is in contrast with that calculated in the paramagnetic (PARA) 
case 0]. A rather strong mass enhancement is retained in the n / 1 case, in 
direct analogy to the paramagnetic case 0. These results are obtained within 
the slave-boson approach in the saddle-point approximation, which we compare 
with the corresponding analysis in the infinite-dimension limit In particular, 
we introduce the concept of a nonlinear staggered molecular field, which shows 
up as the effective (nonlinear-in-magnetization) magnetic gap, evolving at tem- 
perature T = continuously with increasing U from the Slater (Hartree-Fock) 
gap (~ U) into the Hubbard gap. In connection with this evolution, we single 
out the magnetic and Coulomb parts of the localization energy. These particu- 
lar features resolve explicitly the old question about the difference between the 
Slater and Mott-Hubbard insulators in the sense that only the Mott-Hubbard 
gap survives when antiferromagnetism disappears (i.e. above the Neel tem- 
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perature). Most of the results are contained in the mathematical formulation 
established earlier pQ - Here, we only discuss those points in an explicit 
manner. We believe that these points are relevant to the general physics of cor- 
related systems. This reason is also behind publishing, perhaps belatedly, such a 
simple paper. Additionally, since the effective mass is only weakly renormalized 
for the n = 1 case and it never reaches zero for n ^ 1, the role of the quantum 
fluctuations is relatively not as crucial as for the discussion of the paramagnetic 
state. 

The structure of the paper is as follows. In Section II we provide the analytic 
details of the AF solution, which are followed by the numerical analysis (Section 
III) and the discussion in physical terms in Section IV, where we also compare 
the results with those for a doubly degenerate-band case for n = 1. 



2 Antiferromagnetic solution 
2.1 Saddle-point approximation 

To stress the role of the molecular field coming from the electronic correlations, 
we start from the extended Hubbard model, with intersite exchange interactions 
included, i.e. write 

H = X! + U^2n n n H + U ^n^n n + J^Si • Sj - ^N e . (1) 

i&A i e A jeB i£A 

where the labels A and B represent two interpenetrable sublattices, each con- 
taining half (N/2) of available atomic sites. The first term represents single- 
particle hopping of electrons between the sublattices (nearest neighbors), the 
second and the third express the intraatomic interaction of the same magnitude 
on all sites, the fourth includes the Heisenberg exchange between the sublattices, 
and (—fj,N e ) is the reference energy with /i being the chemical potential, and 
N e (< N ) the total number of fermions. The Heisenberg term has been added 
only to provide an illustration of the concept of the exchange field coming from 
the electron correlations (they will add to one another). 

In the mean-field (saddle-point) approximation for the slave bosons, the ro- 
tationally invariant approach of Li et al. |H] and the original Kotliar-Ruckenstein 
1 formulations can be brought to an equivalent form 7 . Explicitly, the six 
bosons e,p^A^,p^A^ and d appearing in the approach fulfill the following five 

constrains in the mean-field approximation 



e 2 1 ~ 2 



p hi +p hi +d =1 > (2) 

where the subscript ( B ) characterizes the sublattices, and fermion operators 
/ko- and /j[ represent the new (quasiparticle) fermion operators appearing in 
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the theory. The effective Hamiltonian in the saddle-point approximation takes 
the form in reciprocal (k) space 



H - E ^ (/ikcriskcr + f Bk(T fAkc 
kcr 



(4) 



1 



E /ika/^ka \ { a A) - M ~ ^<rJzm B 

kcr V 1 

IZ^Bkcr/skcr ( X^ B) - n+-aJzm A 



N 



NUd 2 + —JzmArriB 



where H con contains the constraints appearing in the theory in the form of 
Lagrange multipliers 



N 



(5) 



N 



-A 



(2A) 



(p% + d 2 ) + Af A) ( PAl + d 2 ) 



+ 



-A< 1S > (e 2 
2 1 



N 



xf B) (p| T + rf 2 ) + (Pk + d2 ) 



d 2 +P B ] +PBI 



(IB) 



The summation over (k) comprises the states in the halved Brillouin zone and 
z is the number of nearest neighbors. The Lagrange multipliers A*- 1 " 4 ' and A^ 1S ' 
correspond to the constraint and, since they are spin symmetric we can put 
x (iA) = A (1B) = A (i)_ The spin.dependent Lagrange multipliers \ {2A) and \ {2B) 
represent spin-resolved constraint ©. In the antiferromagnetic (AF) case the 



magnetic moments tub = (p 2 Bt ~ Pw I ) have the opposite sign, i.e. tua 



—tub = m. This means that the constraints A^ obey the relations A^ 2 ' 4 -' = 
A— fj( 2B ) = \^ 2 \ The narrowing factor q renormalizing the bare band energy £k 
assumes the form 



epAa + dp AS 



epBa + dpBa 



V(l -d 2 ~ P A „) (1 - e 2 - ffl V(l - d2 - S (1 - e 2 - ffl 



(G) 



The Hamiltonian can be diagonalized with the help of the Bogolyubov 
transformation [Q. In effect, the free energy functional of the system (per atom) 
takes the form 



F 

N 



k B T 
N 



EM 

kcr 



1 



-/3(E k -/i e 



k B T 
N 



EM 

kcr 



I 



-Pi-Bk-Me/f) 



(7) 
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- Jzm 2 + Udr + fi- 



N 
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where now 



H con = A (1) (e 2 + p 2 + p\ + d 2 - 1) - A (p 2 + p 2 + 2d 2 ) + Am. (8) 



Here the quasiparticle energies are £k = y (? £ k) 2 + A 2 , with the magnetic gap 

((2) (2)\ (2) (2) 

X K A l — AV J = A| — X] , and the effective chemical potential 

M o// = M - A with A = \ (\f ~ A^ 2) ) . The quantity (5 m = \ [\f - a{ 2) ) = 

iA plays a role of the correlation induced molecular field, since it adds to the 
effective Heisenberg field Jzm/2. In the case J = (taken in the numerical 
analysis) (3 m constitutes the entire gap induced by the magnetic ordering (it 
is the magnetic gap). On the whole, the first two terms in J7J provide the 
contribution to the thermodynamics coming from the single particle excitations 
in the magnetic (Slater) subbands caused by AF ordering and having energies 
±£"k- These quasiparticle energies comprise the effective band narrowing or the 
mass renormalization m*/mo = l/g, and the molecular field (3 ml both to be 
determined in a self-consistent manner detailed below. The field (3 m arises from 
the local constraint ©. Thus, one can say that the molecular field is induced 
by the correlations. 

The functional (J7J must be minimized with respect to all appearing Bose 
fields and the chemical potential /i. Effectively, one can reduce J7J) to the form 
with two variables x = 2 A/ (Wq) and d. For the purpose of simplicity and clarity 
of our points we take constant density of states (p(e) — 1/W for — W/2 < e < 
W/2), for which the ground state energy takes the form 

E<; W -q^+x 2 ~ + 9^(1 - nf + x 2 (9) 



N 4 



+ Ud 2 + Qx<? - jmj W, 

where Eq is the ground-state energy, and j = Jz/W. In this expression we have 
already connected m to x via the relation 



■In 



\Jx 2 + (l-n) 2 - (1 - n) 



six 2 + 1 - 1 



(10) 



For the sake of completness, we list also the explicit form of the band narrowing 
factor for n = 1 : 

2 V 



1 



1 - 2r/+ V(l - 2r]) 2 - m 2 . (11 



Note that the variable x has a physical meaning of the ratio of the Slater gap 
to the renormalized band energy. In other words, it provides a relative strength 
of the effective local field with respect to the renormalized kinetic energy. The 
growing ratio A/(Wq) drives the system towards localization induced by the 
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formation of staggered magnetic moments arranged in two sublattices, whereas 
the growing ratio U /W drives the system towards localization independent of 
magnetic ordering. Therefore, the present formulation will allow us to single 
out the contributions coming from the two factors. The magnetic energy is 
thus measured with respect to the band energy (~ Wq) renormalized by the 
Coulomb interaction. 



2.2 Asymptotic limits: U — > oo and U — > 

We now show that the results obtained above provide correctly a mean-field 
solution in the U — > oo limit. For the sake of simplicity consider the half- 
filled band case with J = and for symmetric form of the density of states, 
p(e) — p(—e). We can then write the ground state energy (per site) in the form 



N 



/o 
de p(e)y/(qe) 2 + A 2 + Ud 2 + Am. (12) 
-W/2 

For large U, the gap A is also large. In that limit the energy has the form 

^ = -^<c 2 >+^ 2 + A(m-l), (13) 

where < e 2 >= f_ W i 2 ^ e p{ e ) e2 ■ Minimizing this expression with respect to A 
we obtain that 

m^l-(|) 2 < £ 2 >. (14) 
Minimization with respect to m and d 2 gives the relations 

< 6 > =1 > ( 15 ) 

and 

< e 2 >= 1, (16) 



UA 

where q x = dq/dx. Combining Eas. H14|) - (|16|l we obtain that q — 2(1 — m)q m , 
and A = (q m / 1<p)U ■ Explicitly, from the fact that q ~ 2?7 / (1 — m) and that 
numerically q ss 1, we obtain that d 2 ~ (1 — m)/2, and A w U/2. Thus finally, 
for the featureless form of the density of states we have 

< e 2 > 1 [W^ 2 



"- 1 - 4 — = e{u) ■ < 17 > 

and 

1 - -'^ <-) 



G 



In other words, in the U — ► oo limit the Hubbard model reduces to the Heisen- 
berg model [2] with the Hubbard gap A = U, since in the mean-field approxima- 
tion the ground-state energy is then given by the kinetic exchange contribution 

El 



N UN 



UN ^ UtiBiA " U(N/2)^^ 
' ' de 6 2 p(e) = -2 V, (20) 



with p(e) = J2k^( e ~ e k)- Also, the magnetic gap reduces to the atomic 
value of the Hubbard gap, as 2A — ► U with < e 2 > /U — ► 0. 

The Hartree-Fock (Z7 — * 0) limit is recovered once one notices that the 
method has been constructed to obtain q = 1 in the weak coupling limit 0. 
Under that circumstance Eq.@ reduce to the usual Hartree-Fock form if we 
assume that now x — 2A/W -C 1. This limit was checked out also numerically, 
but the results are not reproduced here. 



3 Numerical analysis 

In Fig.l we have displayed both the effective magnetic (Slater- type) gap 2 A 
and the Mott-Hubbard gap for the paramagnetic phase, both for n = 1 Those 
characteristics are plotted for the ground state. The chemical potential is then 
/u(T = 0) = A. The Mott- Hubbard gap is expressed through the difference in 
the chemical potential in the paramagnetic case (A = 0) for n = 1.001 (the 
upper part) and for n — 0.999 (the lower part) and was discussed earlier [S]. 
For n = 1 the Slater split-band picture appears for arbitrary small U and A 
increases with increasing U/U c , where U c = 8|e|. In the limit U/U c ~ 1 the 
gap is composed of the Slater and the Mott-Hubbard parts, and when — > oo 
the former merges gradually with the latter. This can be seen explicitly in 
Fig. 2, where we have shown the ground state energy E/(WN) versus U/U c . In 
the strong-correlation limit the energy is determined by the kinetic-exchange 
contribution ~ 1/U as shown explicitly. The fitted coefficient a = 0.042 to 
the numerical results cannot be attributed to any particular 3d structure, since 
we have used in numerical calculation a constant value of p(e). The energies of 
para- and antiferro- magnetic states merge for n = 1 and in the U — > oo limit. 
The inset illustrates another interesting characteristic of the solution, namely 
the magnetic gap is not proportional to the magnetization, as one would expect 
from the Hartree-Fock (Slater) solution. In other words, the molecular field (3 m 
is a nonlinear function of m, since from the condition df/dA = we obtain the 
relation 2A/W — — m 2 . This is the reason why the antiferromagnetic 
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U/Uo 



Figure 1: The Slater gap 2 A/14 7 and the Hubbard gap as the functions of 
interaction strength. The Slater gap merges with the Mott-Hubbard gap as 

U/U c -» co. 



u/ u c 




Figure 2: Ground state energy for AF state and n = 1 versus U/U c . The 
lower curve is the fit to the expression (—0.042 x U c /U). The inset displays the 
difference in behavior of magnetic moment m and half of the Slater gap (A/W), 
both plotted as a function of the band filling. 
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2 4 6 8 10 

U/Uc 



Figure 3: The Slater gap parameter A/W (top panel) and the magnetic moment 
m = (riii ~ n ii) (bottom) versus U/U c and for different values of n. 



gap for almost localized fermions cannot be regarded as just the Slater gap. 
Also, the AF solution disappears at the critical band filling n = n c ~ 0.83. 

To visualize the difference between the magnetic gap and the magnetization 
we have plotted in Fig. 3 both quantities as a function of the interaction strength 
U/U c = U/2W for different band filling parameter n. While for n = 1 the 
magnetic moment saturates gradually with growing U/U c , A keeps increasing. 
The magnetic gap for n = 1 increases and eventually A ~ U ; this circumstance 
indicates again that the magnetic gap merges with the Hubbard gap, which can 
be estimated analytically and is « -W + AUW 2 /(z(U 2 + W 2 )) -> U(l - W/U) 
for W/U -> 0. 

The double occupancy probability d 2 = (n^nii) is shown in Fig. 4 for differ- 
ent band fillings. It decreases continuously with growing U/U c , i.e. the charge 
fluctuations are gradually suppressed, while the magnetic moment grows (cf. 
Fig. 3). The difference in the behavior of d 2 and m is caused by the circumstance 
that the d 2 is of intraatomic nature, whereas m is determined from the compe- 
tition between the magnetic energy ~ (3m 2 (also of intraatomic nature) and the 
renormalized band energy ~ Wq . The inset to Fig. 4 exemplifies the difference 
between the diminution of d 2 with growing U/U c for n = 1 in two situations: 
For the paramagnetic (PARA) case d 2 = for U > U c ; this feature is concurrent 
with the well-known effective mass divergence at the Mott-Hubbard localization 
boundary (the Brinkman-Rice point). This divergence does not emerge in the 
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Figure 4: The double occupancy d 2 = (rhfrkl) vs U/U c , and for the different 
n values. The inset display the difference in behavior for para- and antiferro- 
magnetic cases for n = 1. 



antiferromagnetic state as d 2 approaches zero gradually, in the same manner, 

as m approaches saturation m — > 1). 

In Fig. 5 we have displayed the inverse band narrowing factor g _1 (for n = 1) 

which turns into the effective mass renormalization m*/mo = 1/q (for n < 1) 

as a function of /U c and for different n values. Again, the inset illustrates the 

difference with the n = 1 case. One should note that the enhancement factor 

in AF state is very small compared to that in PARA state, which is equal to 
r 2 i -1 

1/q = 1 — (U/U c ) . The difference between AF and PARA states dimin- 
ishes with decreasing n, as in that situation the magnetic moment is reduced 
rapidly. So, the weak band narrowing in the n = 1 case can be associated with 
the circumstance that the Fermi level falls in the gap. This is the reason why 
m* /mo raises rather steeply around n ~ 0.95. Thus, the appearance of the 
itinerant antiferromagnetism with an isotropic gap changes the Brinkman-Ricc 
scenario for the Mott-Hubbard transition, as it has been underlined before ||]. 
Also, the physical parameters d 2 ,l — m and 1 — n are all of the same magnitude. 
This is easy to envisage when estimating e.g. the band narrowing q, which is in 
the AF state roughly sim2d 2 / (1 — m 2 ) and is of the order of unity. 

In Fig. 6 we have displayed the stability regime (n, U/U c ) of the AF phase. 
The full circles has been obtained ^H] m the limit of infinite dimension with 
the help of quantum Monte Carlo simulation. One should note the range of the 
filling n of stable AF phase is the broadest for U/U c ~ 1, i.e. when the molecular 
field is the strongest (cf. Fig. 3). Note also that the Monte-Carlo results did note 
provide the asymptotic behavior for n — > 1, as it does not reduce correctly to 
the Hartree-Fock limit. 
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Figure 6: The stability regime of AF solution; the inset: results of Monte-Carlo 
calculations in the d — >limit [8] 
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Figure 7: Ground state energy vs U/U c for different band filling n. In each case 
the energy is - U for U < O.W c , and ~ 1/U for U > U c . 



For the sake of completeness we display in Fig. 7 the ground state energy 
as a function of U/U c , for different n close to the half filling. In each case (for 
n < 1) the system contains the contribution <~ U for /U c < 0.5 and the contri- 
bution <~ 1/U for U/U c > 1. In other words, the solution interpolates between 
the Hartree-Fock and the kinetic exchange limits. Moreover, the shift of the 
curves with diminishing n in the U — > oo limit is due to the presence of 
the band term proportional to n(l — n). One should note that in the param- 
agnetic state the direct Coulomb interaction contribution (Ud 2 ) competes with 
band energy (—Wq/4). In the antiferromagnetic state the local magnetic con- 
tribution (— |/3 m |m) is of the opposite sign than the Coulomb part. However, 
it changes also the band energy. To specify the role of the staggered field we 
have plotted in Fig. 8 the difference (Ud 2 — \(3 m \m)/W as a function of U/U c . 
The magnetic contribution exceeds the Coulomb part for U > 0.5U C (i.e. away 
from the Hartree-Fock limit), and is particularly strong (and comparable to 
the band energy) in the regime U ~ U c , where the Mott-Hubbard localiza- 
tion would take place for the paramagnetic case. This circumstance indicates 
why the metal-insulator transition present at T = in the paramagnetic phase 
is wiped out when the itinerant AF sets in. Simply, the magnetic contribution 
(<~ /3 m m) freezes gradually the spatially separated electrons in the antiferromag- 
netic phase, taking also advantage of nonvanishing kinetic energy (the kinetic 
exchange, ^< e 2 /U, contribution!). In the PARA phase we have instead a direct 
competition between band and Coulomb energies. 

To summarize, the Mott-localization is achieved gradually at T = in the 
AF state. In other words, the present approach provides a crossover behavior 
from Slater to Mott insulator, as shown e.g. in Fig. 2. The same holds true 
even when we include the intersite exchange (J > 0), which may originate from 
the quantum Gaussian fluctuations. The continuous evolution with growing 
U /Wdoes not preclude the first order transition at nonzero temperature, as has 
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Figure 8: Comparison of the magnetic (— (3 m m) and the Coulomb (Ud 2 ) terms 
in the ground state energy versus U/U c . 



been demonstrated some time ago for the paramagnetic state ^1] and subse- 
quently reconfirmed in the limit d — > oo ^2|- The correlated state builts in 
gradually with increasing U/U c , as can be seen from the gradual increase of the 
mass enhancement in n < 1 case, i.e. when the ground state is always metallic. 
The same holds true for the ground state energy shown in Fig. 7. 

4 Discussion and conclusions 

We have addressed the question of crossover from Slater (Hartree-Fock) to Mott- 
Hubbard (atomic) picture in the half-filled band case, as well as have analyzed 
the behavior of quasiparticle physical properties in AF state in the half- and non 
half-filled-band cases. Although our analysis is based on the saddle-point solu- 
tion within the slave-boson functional-integral approach, the results can serve as 
a mean-field analysis, since they interpolate between those in the Hartree-Fock 
approximation in the limit U — > and those in mean-field approximation for 
the Heisenberg model (for n = 1) in the U — > oo limit. They also represent 
basis for inclusion of Gaussian fluctuations ^Sj m a magnetically ordered (AF) 
state close to the Mott-Hubbard localization. However, one should realize that 
for the half-filled band case the effective mass is only weakly renormalized so 
the renormalization factor Zk F — (mo/m*) ~ 1. Also, for n =/= 1, Zk F remains 
always finite so, perhaps the role of the quantum fluctuations is not as crucial 
for AF state, as it is for the paramagnetic state. In any case, it will be much 
more involved that in the paramagnetic case ^Sj • The full analysis of the Mott- 
Hubbard boundary should also include the disordered local-moment phase |T3] > 
so far not included within the present scheme. 

The physical meaning of the results obtained within the slave-boson ap- 
proach (SBA) can be characterized as follows. It is well known [Q, [E] - [E], 
[T3] that the saddle point solution of the approach reproduces the results of 
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the Gutzwiller approach (GA), as far as the overall (ground-state) behavior is 
concerned. Moreover, SBA leads to an improvement of the GA by incorpo- 
rating both the quasiparticle picture of those systems and the quantum Gaus- 
sian fluctuations (not considered in this paper). In general, the pseudo-fermion 
fields {/io-} representing the quasiparticle states are in one-to-one correspon- 
dence to the original fermion fields in the physical Fock space. The mapping is 
quite obvious within the above analysis. Specifically, the quasiparticle energy 

-Ek = \J (<?£k) 2 + A 2 leads to the following density of states 

2 \E 



p(E) = p°{ 
Hy ' \/E 2 — A 2 V 



(21) 

e=^E 2 -A 2 



where p^(e) is the bare density of states in PM phase. Thus the enhancement 
factor due to the correlations is distinct (and disappear in the first factor) from 
the change of the density of states caused by the appearance of the magnetic 
gap. By analogy, for the paramagnetic state (A — * 0), for which E^ can be 
written as 

E k = qe k = e k + £(e k ), (22) 

and where the self-energy S(w) = (1 — q)uj leads to the mass enhancement 
m*/m = (1 — dY./duj)^ 1 = 1/q, we cannot single out the factor (1/q) in l|21|) as 
the corresponding enhancement also in AF state. In the doubly-degenerate-band 
case (for n = 1) and under the same- type of approximation scheme, the role of 
the magnetic gap is played by the gap formed by an alternant orbital ordering 
in the ferromagnetic state The gap (6) in the latter case is reproduced in 
Fig. 9 for the magnitude of the intraatomic (Hund's-rule) exchange J/U = 0.2. 
Again, the gap approaches asymptotically the Mott-Hubbard gap, which in this 
case [Sj is equal to (U — 3J)/2. However, here there is a critical value of U /W, at 
which the system exhibits antiferromagnetic orbital ordering. In the half-filled 
case for the doubly degenerate case (n — 2) the gap appears for an arbitrary 
small U ^2] in the AF case and for the critical value of U for the paramagnetic 
phase [T5|. 

One should also note that within the present SB scheme the antiferromag- 
netic Slate-type state (AFS) evolves for n = 1 gradually into antiferromagnetic 
Mott (AFI) state with increasing U /W ratio. For n < 1, we observe only antifer- 
romagnetic metallic (AFM) state. This is in contrast with the results obtained 
with the composite-operator method |16| , as well as with an original analysis |17j 
where a phase border line between itinerant- and localized-staggered moment- 
bearing phases is drawn. Localized-moment phase for n = 1 — S, 5 « 1 can 
be obtained only when the polaronic effects due to antiferromagnetic moments 
surrounding the hole in the Mott insulator are included ^SI- Obviously, such 
ASF-AFI border line at temperature T > (for n = 1) is present and is 
induced by the difference in entropy of AFS and AFI states. 

The computational details of the present work are available on request |19| . 
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u/w 

Figure 9: The half-gap parameter S (in units of the bare bandwidth W) vs 
U /W for a quarter-filled doubly-degenerate band. The Hund's rule exchange is 
J = {).2U. 
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Figure Captions 



• Fig.l. The Slater gap 2A/W and the Hubbard gap as the functions of 
interaction strength. The Slater gap merges with the Mott-Hubbard gap 
as U/U c — > oo. 

• Fig. 2. Ground state energy for AF state and n = 1 versus /U c . The lower 
curve is the fit to the expression (—0.042 x U c /U). The inset displays the 
difference in behavior of magnetic moment m and half of the Slater gap 
(A/W), both plotted as a function of the band filling. 

• Fig. 3. The Slater gap parameter A/W (top panel) and the magnetic mo- 
ment m — (njf — riii) (bottom) versus U/U c and for different values of 
n. 

• Fig. 4. The double occupancy d = (n^riii) vs U/U c , and for the different n 
values. The inset display the difference in behavior for para- and antiferro- 
magnetic cases for n = 1. 

• Fig. 5. The effective mass enhancement to* /too (with respectto the band 
value Too) vs U/U c and for the n values shown. The inset shows a rather 
weak enhancement close to the Mott-Hubbard limit. 

• Fig. 6. The stability regime of AF solution; the inset: results ofMonte- 
Carlo calculations in thed — > limit 

• Fig. 7. Ground state energy vs U/U c for different band filling n. In each 
case the energy is ~ U for U < 0.5U C , and ~ 1/U for U > U c . 

• Fig. 8. Comparison of the magnetic (—f3 m m) and the Coulomb (Ud 2 ) terms 
in the ground state energy versus U/U c . 

• Fig. 9. The half-gap parameter 5 (in units of the bare bandwidth W) 
vs U/W for a quarter-filled doubly-degenerate band. The Hund's rule 
exchange is J = 0.2U. 
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